arXiv:1507.02333v2 [astro-ph.SR] 29 Oct 2016 


DRAFT November 1, 2016 

Preprint typeset using DTTjK style emulateapj v. 5/2/11 


LINKING STELLAR CORONAL ACTIVITY AND ROTATION AT 500 MYR: A DEEP CHANDRA 

OBSERVATION OF M37 

Alejandro Nunez\ Marcel A. Agueros\ Kevin R. Covey^, Joel D. Hartman^, Adam L. Kraus*, 

Emily C. Bowsher*, Stephanie T. Douglas*, Mercedes Lopez-Morales®, David A. Pooley®’'", Bettina Posselt®, 

Steven H. Saar®, Andrew A. West® 

DRAFT November 1, 2016 

ABSTRACT 

Empirical calibrations of the stellar age-rotation-activity relation (ARAR) rely on observations 
of the co-eval populations of stars in open clusters. We used the Chandra X-ray Observatory to 
study M37, a 500-Myr-old open cluster that has been extensively surveyed for rotation periods (Prot)- 
M37 was observed almost continuously for five days, for a total of 440.5 ksec, to measure stellar 
X-ray luminosities (Lx), a proxy for coronal activity, across a wide range of masses. The cluster’s 
membership catalog was revisited to calculate updated membership probabilities from photometric 
data and each star’s distance to the cluster center. The result is a comprehensive sample of 1699 
M37 members: 426 with Prot, 278 with X-ray detections, and 76 with both. We calculate Rossby 
numbers, Ro = Prot/^, where r is the convective turnover time, and ratios of the X-ray-to-bolometric 
luminosity, Px/Pboi, to minimize mass dependencies in our characterization of the rotation-coronal 
activity relation at 500 Myr. We find that fast rotators, for which Ro < 0.09 ± 0.01, show saturated 
levels of activity, with log(Px/Pboi)= —3.06 ± 0.04. For Ro > 0.09 ± 0.01, activity is unsaturated 
and follows a power law of the form Pf, where /3=— 2 . 03 ^ 0 / 4 . This is the largest sample available 
for analyzing the dependence of coronal emission on rotation for a single-aged population, covering 
stellar masses in the range 0.4-1.3 Mq, Prot in the range 0.4-12.8 d, and Lx in the range io28-4-30.5 
erg s“*. Our results make M37 a new benchmark open cluster for calibrating the ARAR at ages of 
«500 Myr. 

Subject headings: Galaxy: open clusters and associations: individual (M37) - stars: activity - stars: 
coronae - stars: rotation - X-rays: individual (M37) 


1. INTRODUCTION 

X-ray emission traces the strength of a stellar magnetic 
dynamo and is tightly linked to stellar rotation and age. 
Rotational shear at the boundary between the radiative 
and convec tive zones is thought to power the solar dy¬ 
namo (e.g., Parker|[l993 l. The generated magnetic fields 
heat the corona to temperatures sufficient to produce X 
rays (|Parker|1979| |Cowling|1981| . Observations of solar¬ 
like main-sequence (iVib) stars confirm this connection 
between rotation and X-ray emission, finding that the 
X-ray luminosity Lx fx v^ot^, where Vmt is the stellar 
surface rotational velocity and /3 ~ 2 (Pallavicini et al. 
1981 Pizzolato et al.|[2003 l. 


Tn addition, observations have shown that stars rotate 
more slowly as they age, with rotation decaying with 
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time such that Urnt, (x t with a = 0.5 for solar -mass 


stars (e.g., Skumanich 1972 
gular momentum lost throug 


Meibom et al. 20111. An- 
1 stellar winds is generally 


thought to be responsible for the Skumanich law, with 
the exact dependence of v^ot on age being a function of 
the assur ned field geometry and degree of core-enve lope 
coupling ( Kawaler| 1988 Krishnamurthi et al.||1997). 

Combined, these two effects predict that for solar- 
type stars. Lx should decay such that Lx oc with 

a/3 = 1. It is also known, however, that this picture of 
a uniformly decaying Lx does not apply to stars with 
ages ^100 Myr. Surveys of solar-type stars in the Orion 
Nebula and the Pleiades (^10 and 100 Myr, respectively) 
found that Lx decays more slowly: Lx ex ( Queloz| 

et al. II 1998 Preibisch fc Feigelson||2005 1. Because of the 


paucity of older, nearby clusters, constraints are more 
limited at ^200 Myr and come mostly from field stars 
(e.g., Gudel|[2007 1 : observations of five sol ar analogs in¬ 
dicat ed that Lx oc for t > 1 Gyr (Giidel et aLf 

19971), as did a survey of the 1-2 Gyr cluster INCC 752 


(jciardino et al. 2008). A change in the magnetic field 
topology or core-envelope decoupling are the commonly 
invoked explanations for this sharp drop-off in Lx, but 
it remains poorly understood. Chromospheric activity, 
also thought to trace the strength of the stellar dynamo, 
may also suffer a simila rly steep decline at ages >1 Gyr 
dPace & Pasquini||2004|). 

T he homogeneous, co-eval populations of open clusters 
provide an ideal environment for calibrating empirically 
the stellar age-rotation-activity relationship (ARAR). 
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The Hyades’s proximity («50 pc) has facilitated obser¬ 
vations of its members at many wavelengths and estab¬ 
lished it as the benchmark mi ddle-aged (500-600 Myr) 


cluster (Perryman et al. 1998). Empirical calibrations 
of the AKAK often interpolate directly between the ob¬ 
served properties of solar-type Hyads and that of the Sun. 
The Hyades therefore serves as the only anchor for our 
understanding of coronal emission between 500-600 Myr 
and 4.5 Gyr, a significant fraction of a solar-type star’s 
MS lifetime. 

However, the Hyades’s proximity is a drawback for X- 
ray studies. The cluster spans well over 10°, and mod¬ 
ern X-ray telescopes with small fields of view are poorly 
suited to large-scale studies of the cluster. Indeed, the 
most recent such study was conducted by [Stern et al.| 
(|1995|), who detected 185 Hyads with ROSAT. 


Recently, Douglas et al. (2014) probed the rotation- 
activity relation tor stars in tne Hyades and Praesepe, an¬ 
other «600-Myr-old cluster, using both coronal and chro¬ 
mospheric (Thq) activity measurements. These authors 
found that coronal activity declines faster than chromo¬ 
spheric activity as a function of (unsaturated) Rossby 
number = Prot/'^j where r is the convective turnover 
time. However, because of the paucity of stars in either 
cluster with measured rotation periods (Prot) and X-ray 
detection (24 in the Hyades and 15 in Praesepe), the ex¬ 
act relationship between coronal activity and rotation at 


Other studies (e.g.. 

Pallavicini et al.[ 1981| 

Giidel 

et al.|1997||Randich|2000[[iVlamaiek & Hillenbrand 

rpnnsf 

Wright et al. 201ip have examined the stellar rotation- 


lar issues: either there are too few stars at a given age to 
get significant statistics or the population is too hetero¬ 
geneous to strongly constrain age-dependent properties. 

M37 (NGC 2099), a Hyades analog at Ril.5 kpc, has 
been extens ively surveyed in the optical by Messina] 


et al.j (|2008|) (hereafter ME08), [Hartman et al.j (|2008a 
(hereafter HA08), and Hartman et al. 1 2009 r(hereafter 


HA09), resulting in hundreds of Ppot measurements. As 
part of our Palomar Tran sient Factory Open Gluster Sur¬ 
vey (Agiieros et al. ]MT] ), we obtained a complementary 
deep (440.5 ksec) Chandra observation of the cluster. 
M37 is distant enough that a significant number of clus¬ 
ter stars fall within the footprint of Chandra’s AGIS-I 
chips. Our X-ray data, combined with M37’s rich mem¬ 
bership and the extensive set of available P^t data, make 
it an excellent laboratory for examining the dependence 
of coronal activity on rotation. 

We begin in Section by describing our Chandra data 
and data reduction, and the construction of our X-ray 
catalog. In Section we consolidate the ME08 and 
HA09 Prot data for the cluster, and we re-visit the cluster 
membership using our own membership probability cal¬ 
culations (fully described in Appendix |A]) In Section]^ 
we derive stellar quantities using optical^nd X-ray data, 
and we match our X-ray sources to the consolidated list 
of optical objects. In Section we present our results 
and compare them to those in the literature. We con¬ 
clude in Section We present spectral and timing anal¬ 
ysis of our highest-count X-ray Chandra sources in Ap¬ 
pendix 1^ 

2. CHANDRA OBSERVATIONS 



Fig. 1. — Thick blue squares (with chip numbers) indicate the 
footprint of the Chandra ACIS chips for four of our observations; 
the thin orange squares indicate that of our fifth observa t ion. A 
40^x40^ Megacam image of M37 from |Hartman et al.| (|2008a|| 
serves as the background. 


TABLE 1 

Log of Chandra Observations 


ObsID 

MJD 

Exposure^ 

Nominal Aimpoint^ 

Roll 



(s) 

“J2000 

<^J2000 

(°) 

13655 

55879 

51179 

05:52:17.61 

-1-32:33:59.0 

103 

13656 

55886 

165493 

05:52:17.61 

-132:33:59.0 

103 

13657 

55941 

92172 

05:52:18.82 

-132:33:05.0 

253 

14366 

55882 

65074 

05:52:17.61 

-132:33:59.0 

103 

14367 

55883 

60813 

05:52:17.61 

-132:33:59.0 

103 


^ Exposure times are the net usable times after various filtering 
steps are applied in the data reduction process. 

^ The time-averaged location of the optical axis. The aimpoints 
are obtained from the satellite aspect solution before astrometric 
correction is applied. Units of right ascension are hours, minutes, 
and seconds; units of declination are degrees, arcminutes, and arc- 
seconds. 


The central field of M37 was observed five separate 
times between 2011 Nov 14 20:58 and 2011 Nov 19 15:31 
UTG for a total of 440.5 ksec with the Advanced GGD 
Imaging Spectrometer (ACIS; Garmire 2003). The four 
ACIS-I chips and the ACIS-S3 chip were used in Very 
Faint telemetry mode to improve the screening of back¬ 
ground events a nd thus increas e the sensitivity of ACIS 
to faint sources (Vikhlinin 2001). The exposure-weighted 
average aimpoint of the 16:9xl6f9 ACIS-I field-of-view is 
a = 05‘'52“17."86, (5 = -h32°33'48"23 (J2000). The pitch 
angle for four observations was 103°; due to scheduling 
constraints, it was 253° for the fifth. Table provides 
the basic information for our five Chandra observations 
and Figure [2 shows their footprints superimposed on a 
40'x40' i' image centered on M37 obtained by HA08 with 
the Megacam on the MMT telescope. 

The data reduction was done using th e Chandra Inter- 


active Analysis of Observations (CIAO; Fruscione et al. 
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2006 we used CIAO 4.6 and CALDB 4.6.1.1) tools. We 
began by using chandra_repro, which executes all rec¬ 
ommended data processing steps on the level 1 data, in¬ 
cluding subpixel randomization using the EDSER algo¬ 
rithm. This tool creates a new bad pixel file by identi¬ 
fying known bad pixels and pixels with bad bias values, 
searching for afterglows and hot pixels, and marking pix¬ 
els adjacent to afterglows and hot pixels. chandra_repro 
also creates a new level 2 event file that includes only 
events with grades 0, 2, 3, 4, 6, and a status of 0. 
To avoid excluding potentially good events in modestly 
bright point sources, we did not use the very-faint cor¬ 
rection in chandra_repro. 

2 .1. X-ray Source Extraction 

Our source extraction procedure was divided into two 
steps. First, we merged our observatio ns using CIAO and 
ran wavdetect (Freeman et al.||2002|) with a liberal sig¬ 
nificance threshold to obtain a raw candidate source list. 
Second, we used the ACIS Extrac t (AE) point-source 
analysis software (Broos et al. ]Mo| we used AE Version 
2014febl7)P^to prune spurious detections from our can¬ 
didate list and extract photometry. AE’s algorithms are 
tailored to calculate source properties accurately when 
multiple observations are merged and analyzed. 


2 . 1 . 1 . Merging the Observations 

To construct the merged event file, we first ran 
wavdetect on the observations using a false-positive 
probability threshold of 10“®. We then used 
reproject_aspect to register the observations to the as¬ 
trometric frame of the longest exposure (Obsid 13656), 
using only source regions with SRC_SIGNIFICANCE > 7, 
a 3" matching radius, and a residual rejection limit (a 
parameter used to remove source pairs based on pair 
positional offsets) of 076. The values of these parame¬ 
ters were chosen to have reprojec t aspect work w ith a 
more stringent source list (see e.g.,jXue et al.|2011|). The 
linear translations ranged from —Cr.'38 to 0763, rotations 
ranged from —0°.006 to 0°.018, and scale changes ranged 
from 0.9997 to 1.0002. 

We used the CIAO tool merge _obs to produce a 
merged event file, a combined exposure map, and 
exposure-corrected flux images. We constructed flux im¬ 
ages and exposure maps for three bands: 0.5-7 keV (full), 
0.5-2 keV (soft), and 2-7 keV (hard). 

Figure [^shows the full band effective-exposure map for 
the mergra 440.5 ksec event. Since one of our five obser¬ 
vations is rotated 150° relative to the others, the merged 
exposure map is inhomogeneous, with the longest ex¬ 
posure near and around the average aimpoint. We used 
ds9 to produce a representative color composite image by 
combining, overbinning, and smoothing (using a Gaus¬ 
sian of kernel radius 5) the 0.5-1.5 keV (red), 1.5-2.5 keV 
(green), and 2.5-7 keV (blue) exposure-corrected flux im¬ 
ages, and the result is shown in Figure 

2 . 1 . 2 . Source Detection and Pruning 

To obtain an initial list of candidate X-ray sources, we 
first created point-spread function (PSF) maps for the 

Available from http://www2.astro.psu.edu/xray/docs/ 
TARA/ae_users_guide.html 



88.4 88.3 88.2 88.1 88.0 87.9 87.8 87.7 


RA (J2000) 

Fig. 2.— Effective exposure map for the full (0.5-7 keV) band of 
our merged Chandra event displayed in logarithmic gray scale (in 
units of ksec). The average aimpoint is indicated by the plus sign. 


individual observations using mkpsfmap with an encircled 
counts fraction (the fraction of the PSF to enclose) of 0.9. 
We then selected the minimum PSF map size out of the 
five observations at each pixel using dmimgfilt, which 
allowed us to detect point sources that are smaller than 
the mean size, but still larger than the local PSF in the 
individual maps. 

To obtain an initial candidate source list, we ran 
wavdetect using this minimum PSF map and the merged 
event file and exposure map. We used a scale sequence 
ranging from 1 to 8\/2 for the ACIS-I chips and from 8 to 
P^lin both cases the scal es had factor 
(20011 pointed out 


Alexander et al. 


32^2 for ACIS-S3 
increments of \/2. 
that real sources can be inadvertently lost using a too 
stringent false-positive probability threshold. We there¬ 
fore reset our threshold to 10“® in wavdetect to allow 
for a non-negligible number of spurious sources to be in¬ 
cluded in the initial candidate list. 

Separate candidate source lists were prepared for the 
full, soft, and hard energy bands. To consolidate the 
three lists, we matched these lists against each other us¬ 
ing a 2" matching radius for sources with angular sep¬ 
aration 9 < b' from the average aimpoint, and 4" for 
sources with 9 > b'. No multiple matches occurred us¬ 
ing this approach. Once matched, the central location of 
each point source was adopted, in order of priority, from 
the full, soft, or hard energy band list. The resulting list 
consists of 794 candidate X-ray sources: 694 detected in 
the soft band (53 only in the soft band), 173 in the hard 
band (10 only in the hard band), and 608 in the full band 
(50 only in the full band). 

To prune spurious sources from this list, we used AE’s 
binomial probability P (?, Appendix A2)j^ P is the 
probability that, based on the source anoTocal back¬ 
ground measurements, a candidate source is a back¬ 
ground fluctuation. This approach to producing a final 
source catalog is preferable to implementing wavdetect 
with a more stringent false-positive probab i lity threshold 
(e.g., Lehmer et al.||2009 Xue et aLl|2011 Wright et al. 


2014 li Not only is the pruning based on a transpar- 


ent mathematical criterion, P, but AE also uses a more 
complex approach to treat source regions in multiple im- 


Larger scales are needed for ACIS-S3 because of the large 
PSFs far from the aimpoint. 

prob_no_source, see section 5.10.3 of the AE manual, http: 
//www2.astro.psu.edu/xray/docs/TARA/ae_users_guide.pdf 
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Fig. 3. — Chandra smoothed RGB composite image of the central region of M37. We combined, overbinned, and smoothed (using a 
Gaussian of kernel radius 5) the 0.5-1.5 keV (red), 1.5-2.5 keV (green), and 2.5—7 keV (blue) exposure-corrected flux images from our 
merged Chandra event file. The plus sign indicates the average aimpoint. 


ages by replacing wavdetect’s ellipses with polygons that 
more closely simulate the PSF, and by calculating more 
accurate source positions, which maximizes the signal- 
to-noise ratio. We followed AE’s suggestion and rejected 
candidate sources with P > 0.01 as spurious. 

To refine the source positions, we improved those cal¬ 
culated by wavdetect using AE’s source positions esti¬ 
mates. For source regions with 0 < 5' we used the “mean 
data” position estimate, which is simply the mean posi¬ 
tion of events. For source regions with 0 > 5', we used 
the “correlation” position estimate, which determines the 
source center by correlating the neighborhood around the 
source (not just the extracted counts) with the PSF of 
the source. This is a more accurate source center posi¬ 
tion when the PSF is asymmetric, which is the case for 
source regions far off-axis 

The absolute astrometry of the individual observations 
was then corrected by using reproject_aspect again to 
register them to the astrometric frame of the 2MASS 
Point Source Catalog ( Cutri et al. [2003 1. This time, the 
procedure used the refined source positions of our X-ray 
sources and a 2MASS list of point sources of photometric 
quality in the JHKs bands within a square of side 
20' centers on the ACIS-I average aimpoint («1900 ob¬ 
jects). We used a 2'.'5 matching radius and a residual 
rejection limit of C/.'fi. When multiple matches occurred, 
the closest one was selected as the true match. The linear 


See section 5.3 of the AE manual, http://www2. astro. 
psu.edu/xray/docs/TARA/ae_users_guide.pdf tor a discussion of 
source center estimates. 

Source confusion was not an issue for the 2MASS sources in 
the list we used. 


translations ranged from —0'.'28 to —O'.'16, and a rotation 
angle of —0°.01 and scale factor of 1.0003 was applied to 
all five observations. 

We iterated over the process of source and local back¬ 
ground extraction, P calculation, candidate source list 
purging of P > 0.01 sources, source position refine¬ 
ment, and absolute astrometric correction, until no more 
sources were discarded. This was necessary because 
source and local background extraction in surviving 
sources may be affected by the disappearance of nearby 
purged sources. Our final Chandra catalog contains 774 
sources. Of these, 321 were detected in the three en¬ 
ergy bands, four were detected in the full band only, and 
twelve in the soft band only. Almost half of the discarded 
sources were originally detected only in the full band. 

The regions used to extract these sources are shown 
in Figure drawn on the merged full-band flux image, 
smoothed using a Gaussian function of kernel radius 4. 
Figure 1^ shows the distribution of the binomial proba¬ 
bility for all the sources in our catalog. This histogram 
clearly shows that most sources have negligible proba¬ 
bilities of being false detections. The final merged event 
file has an absolute astrometry accuracy of <0'.'63, based 
on 210 matches to 2MASS. We investigate the positional 
uncertainty of the individual X-ray sources in more detail 
in Section 12721 

2 . 1 . 3 . Extracting X-ray Photometry 

We followed AE’s recipe for photometry extractiorpq 
using the three energy bands (full, soft, and hard) d^ 

See http://www2.astro.psu.edu/xray/docs/TARA/ae_users_ 
guide/procedures/photometry_procedure.txt 
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Fig. 4.— Merged full-band (0.5-7 keV) flux image smoothed using a Gaussian function of kernel radius 4, with the polygons used for 
extracting source events, photometry, and spectra drawn in blue. The average aimpoint is indicated by the plus sign. Red stars indicate 
the four bright sources for which we perform spectral and timing analysis in Appendix]^ We identified 774 X-ray sources. The varying 
sizes of the source regions illustrate the degrading of the PSF far from the aimpoint. 



log P 


Fig. 5. — AE’s binomial probability P that a source is just a 
background fluctuation for all sources that have P < 0.01 in at 
least one of the three energy bands. We considered sources with 
P > 0.01 to be spurious detections. For clarity, P values lower 
than 10“^*^ are set to 10“^*^. 

scribed in Section 12.1.11 AE first defines each extrac¬ 
tion region with a polygon that approximates the «90% 
encircled-energy fraction contour of the local PSF mea¬ 
sured at 1.497 keV. It then calculates exposure times, 
source counts, and background counts for each region 
and energy bandf^ 

AE executes these steps on each Chandra observation 
independently and then merges the extraction data from 
each source by optimizing the source’s signal-to-noise ra¬ 
tio. AE then computes photon fluxes using net (i.e.. 

For background counts, we used AE’s BETTER_BACKGR0UNDS 
algorithm, which models several background components simulta¬ 
neously. See section 7.6.1 of the AE manual, http://www2. astro. 
psu.edu/xray/docs/TARA/ae_users_guide.pdf 


IZZIO.5-7 keV 



1 3 10 30 100 300 


Net counts 

Fig. 6. — Cumulative distribution of sources in our catalog based 
on their net counts in the full (0.5-7 keV), soft (0.5-2 keV), and 
hard (2-7 keV) bands. 

background-subtracted) source counts, exposure times, 
and ancillary response files generated by the CIAO tool 
mkarf. 

Figures!^ and 1^ show the net counts and absorbed X- 
ray flux clistributions in the three energy bands for all 
the sources in our catalog. We calculated the latter by 
multiplying the band photon flux by its median photon 
energy. The median number of net counts for the full, 
soft, and hard band is ~19, 14, and 3, respectively. There 
are 78 sources with more than 100 net counts in the full 
band. In general, hard band detections are less common 
than full band detections, and most hard band detections 
have fewer than ten net counts. The median absorbed 
X-ray flux for the full, soft, and hard band is 4.8 x 10“^®, 
2.4 X 10“^®, and 3.0 x 10“^® erg cm“^ s“^, respectively. 
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Fig. 7.— Cumulative distribution of sources in our catalog based 
on their absorbed X-ray fluxes in the full (0.5-7 keV), soft (0.5-2 
keV), and hard (2-7 keV) bands. 


2 .2. Estimating the Source Positional Uncertainties 

To parametrize the dependence of the X-ray positional 
uncertain ty on 9 and sour ce counts, we applied the ap¬ 
proach of Kim et al. (20071, who described the positional 
uncertainty as 


logAa; = ao -I- ai6l -I- a2log C, 


( 1 ) 


where is in arcseconds, 0 is the off-axis angle in 
arcminutes, C is the net full-band counts, and oq, oi, 
and 02 are constants. This parametrization takes into 
account the fact that the statistical uncertainty of the 
calculated PSF center of a Chandra source increases 
with lower source counts, and that the size of the 
PSF increases with larger 9. We performed our own 
pa rametrization of th e relation above because the one 
by Kim et al. (20071 was based on wavdetect-derived 
positions and photornetry. 

We started by doing a preliminary match of our X-ray 
sources to objects in the optical survey of HA08, which 
is registered to the 2MASS astrometric frame, using 2'.'5 
and 470 matching radii for X-ray sources with 9 < U 
and > 5', respectively. We selected these separations 
by inspecting a histogram of the number of matches as 
a function of angular separation, and selecting the val¬ 
ues that resulted in a mismatch probability <1%. We 
excluded likely false matches based on full-band counts, 
9, and matching radius (e.g., a source with >100 net 
counts, 9 < 5, and matching offset >1"). When multiple 
matches occurred, the closest one was selected as the true 
one. Using these criteria, we identified 452 X-ra y sources 
as having optical counterparts. In Section |4.3| we use a 
more rigorous matching technique that considers the de¬ 
rived positional uncertainty of each X-ray source from 
this section to create our final list of X-ray-to-optical 
matches. 

Next, we performed a multivariate minimization of 
Equation using the offsets between our X-ray sources 
and HA08 optical objects as proxy for positional uncer¬ 
tainties. Once we obtained values for oq, oi, and 02 , we 
adjusted the value of oq until >80 % of our X-ray sources 
with optical counterparts had A^, values larger than 
the match offsets. The final values are Og = —0.2928, 
oi = 0.1200, and 02 = 0.3401. 

We then used these values in Equation to calculate 
the «80% confidence positional uncertainties for all our 
sources, and we report these values in Table 1^ Figure |8] 



Fig. 8.— Positional offset for our X-ray sources with optical 
counterparts in the HA08 survey, shown as a function of distance 
6 from the Chandra aimpoint. Black stars are X-ray sources with 
>100 net counts in the full band; green squares, between 30 and 
100, orange circles, between 10 and 30; and small gray circles, <10. 
We indicate Ri80% positional uncertainty curves for sources with 
10, 30, and 100 net counts following Equationj^ 


shows the positional offset as a function of 9 for the X- 
ray-to-optical matches just described. We overplot three 
curves generated from Equationusing C = 10, 30, and 
100 counts. 

To estimate the expected number of false matches in 
our catalog, we shifted the X-ray source positions by 30" 
out to 6' in all directions, and then re-matched them 
to the optical positions using a 21'5 matching radius. 
The 2'.'5 value was selected by inspecting a histogram 
of the number of matches as a function of matching ra¬ 
dius (from 0'.'2 to 4'.'0) at all shifted locations. We found 
that the median number of false matches is 96 («12.5% 
of our matched catalog) with a median offset of 1'.'7. 

Table|2]describes the 31 columns in our Chandra source 
catalog. Tablej^is a sample of the full catalog. This table 
is available online in its entirety. 

3. OPTfCAL OBSERVATfONS 
3.1. Description of the Optical Surveys 

ME08 selected images from the survey taken with the 
1.0-m telescope at the Mt. Lemmon Optical Astron¬ 
omy Observatory, AZ, of a 2272x2272 area centered on 
M37, and used both the Scargle-Press and CLEAN peri- 
odogram techniques to measure variability for stars with 
13 < U < 20 mag. ME08 found 135 variables and re¬ 
ported Prot for 120 of them, and the adopted P^t for 
each star is the one from the two periodogram calcula¬ 
tions with the least-scattered phased light curve. The 
typical Prot uncertainty is <2%. Based on (P — V) col¬ 
ors, these authors identified 28 of the objects as K stars, 
60 as G stars, and 32 as F stars. For the latter, ME08 
noted that their periods indicate pulsation instead of ro¬ 
tation, as F stars are not expected to display magnetic 
activity and, thus, no star spots. 
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TABLE 2 

Overview of Columns in the Chandra Source Catalog 


Column 


Description 


1 X-ray catalog sequence number (XID), sorted by RA. XID > 1000 are for sources in the ACIS-S3 chip. 

2 lAU designation. 

3, 4 Right ascension and declination of X-ray source (in decimal degrees) for epoch J2000. 

5 ~80% confidence level X-ray positional uncertainty. 

6 Off-axis angle 6 of X-ray source (i.e., distance to average aimpoint of Chandra observations). 

7-15 Source net counts in the 0.5-7 keV, 0.5-2 keV, and 2-7 keV band and corresponding upper & lower errors. 

16-18 Absorbed energy fiux (full, soft, and hard bands.) 

19 Median photon energy (full band), corrected for background. 

20 Fraction of the PSF (at 1.497 keV) enclosed within extraction region. PSF fraction significantly below 90% may 
indicate a source in a crowded region. 

21 Absorption-corrected 0.1-2.4 keV X-ray luminosity —for sources with optical cluster member counterpart only. 

22 Standard deviation of 0.1-2.4 keV X-ray luminosity —for sources with optical cluster member counterpart only. 

23 Logarithmic probability P that extracted X-ray counts are solely from background. 

24 Photometric significance computed as net counts divided by the upper error on net counts. 

25 Source anomalies: (f) fractional time that source was on a detector is <0.9; (a) photometry and spectrum 
may contain > 10% afterglow events. 

26 Variability characterization based on K-S statistic (full band) from the single Obsid showing the most variability: 

(a) no evidence for variability (0.05 <Pks)) (t>) possibly variable (0.005<PKS<0-05)j (c) definitely variable (Pks< 0-005). 
No value is reported for sources with fewer than four net counts in the full band. 

27 ID of optical counterpart. 

28, 29 Right ascension and declination of the optical counterpart. 

30 Offset between X-ray source and optical counterpart. 

31 Notes on the source. 


TABLE 3 

Chandra CATALOG: BASIC SOURCE PROPERTIES 



Source 


Position 




Extraction 


Properties 

XID 

CXOU J 

a (J2000) 
(°) 

S (J2000) 
(°) 

Error 

(") 

d 

(') 

^f,net 

(counts) 

C^s,net 

(counts) 

C*h,net 

(counts) 

PSF Frac. 

Anom. 

Var. 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(10) 

(13) 

(20) 

(25) 

(26) 

1 

055132.87-1-323928.4 

87.886986 

32.657895 

5.5 

10.9 

25.3 

18.6 

6.6 

0.89 

f. 

b 

2 

055133.77-1-322934.0 

87.890743 

32.492781 

3.4 

10.1 

53.0 

21.4 

31.6 

0.90 


a 

3 

055134.32+324019.9 

87.893001 

32.672222 

4.2 

11.1 

66.8 

51.1 

15.7 

0.91 


c 

4 

055136.34+323835.8 

87.901458 

32.643302 

1.9 

9.8 

247.8 

52.6 

195.2 

0.90 


b 

5 

055136.98+324011.7 

87.904097 

32.669942 

3.8 

10.6 

61.2 

32.9 

28.3 

0.91 

f. 

a 


Note. — This table is available in its entirety in the electronic edition of the ApJ. The first five rows with some data columns are shown 
here for guidance regarding its form and content. Table [^describes all the columns in this table. 


HA08 used g'r'z' filters on Megacam and the 6.5-m 
MMT telescope at the MMT Observatory, AZ, to sur¬ 
vey a 24' X 24' area around the center of M37. They ob¬ 
tained photometry for >38,000 objects in at least one 
band. HA08 converted the photometry of objects that 
had measurements in all three bands (~ 16,500) to gri 
magnitudes. The typical one-sigma error is 0.01 mag 
for the gri HA08 photometry, which covers the range 
15 ;< r ;< 23 mag. 


Hartman et al. (2008b I identified 1445 variable stars in 
the HA08 survey. These authors used the Lomb-Scargle, 
box-fitting least-squares, and phase-binning Analysis of 
Variance (AoV) algorithms to measure periodicity. The 
most likely period for each star was selected from these 
three measurements by eye. 

Finally, HA 09 identified periodic va riability due to ro¬ 
tation for 575 Hartman et al. (2008b I variable stars. For 
these stars, HA09 deter mined Tj-nt. using the multihar - 
monic AoV algorithm of Schwarzenberg-Czerny (19961, 
which fits the harmonic series 


to light curves, where are the amplitudes and tpi are 
the phases. HA09 calculated three Prot using V = 1, 2, 
and 3. For the 372 stars with <10% differences between 
the three U-ot values, the N=2 P^ot was adopted; these 
are referred to as clean peri ods. For the other stars , the 
(non-clean) period quoted in Hartman et al. (2008b I was 
kept. Here the typical P^ot uncertainty is <1^ 


3.2. Consolidating the Optical Surveys 

We consolidated photometry and Prot from the surveys 
described above. We matched the ME08 and HA08 po¬ 
sitions using a 10" matching radius and inspected by eye 
multiple-match cases to determine the most likely match 
using available gri and BV photometry as guidance. The 
liberal matching and individual inspection was necessary 
because, as noted by HA09, the absolute astrometry of 
ME08 objects displays a widespread non-linear distor¬ 
tion. Of the 153 stars in the full ME08 catalog, 16 have 
no HA08 counterpart. We found 37 matched stars to 
have ME08 and HA09 Prot measurements |F For consis- 


N 

f{t) = 00-1-^ aiCOs(27rit/Prot + 4'^), 

i=l 


HA09 listed 22 matches to ME08 with clean Prot, but we find 
that their match between HA09 ID 424 and ME08 ID 3208 appears 
to be a mis-identification based on their available photometry. 
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TABLE 4 

Overview of Columns in Catalog of Optical Objects 


Col. Description 


I Source ID from HA08(>10000), ME08 (<10000), 
or literature. 

2, 3 Right ascension and declination of object (J2000). 

4-9 gri magnitudes and their standard deviations. 

10 Rotation period Prot adopted for object. 

II Rotation period flag: (0) non-clean HA09 period; 

(1) clean HA09 period; (2) ME08 period. 

12-14 Membership probabilities (Pg, Pb, and Pmem)- 

15 Stellar mass. 

16 Convective turnover time r. 

17 Bolometric luminosity Lboi- 

18 Standard deviation of Lboi- 

19 X-ray sequence number (XID) of X-ray counterpart. 

20 Offset between optical object and X-ray counterpart. 

21 ID of ME08 counterpart —for HA08 objects only. 

22, 23 Right ascension and declination of ME08 counterpart. 
—for HA08 objects only. 

24 Offset between HA08 object and ME08 counterpart. 
—for HA08 objects only. 

25 Notes on the object. 


tency, we adopted the HA09 period for all 37. 

We used gri photometry for our analysis. For ME08 
stars, we adopted the gri magnitudes of their HA08 coun¬ 
terparts. Of the 16 ME08 stars with no HA08 counter¬ 
part, we found three with available SDSS photometry. 
For the remaining 13 ME08 stars, we derived gri pho¬ 
tometry from the V and (i3 — V) reported by ME08 u sing 
the transformations derived by Jester et al. (20051 and 
Lupton (2005) ^ 

Our consolidated optical catalog includes 16,584 ob¬ 
jects in the M37 field with gri photometry. Of these, 
657 stars have a Prot measurement. Table summarizes 
the 25 columns in our catalog of optical objects, a pre¬ 
view of which is presented in Table 1^ This table, which 
includes all objects from our optical catalog with either 
a Prot, a non-zero probability of being a cluster mem¬ 
ber (se e dis cussion below), or an X-ray counterpart (see 
Section 4.31, is available online in its entirety. 


3.3. Membership Determination 

ME08 and HA08 relied solely on distance from a fidu¬ 
cial isochrone in a color-magnitude diagram (CMD) to 
determine M37’s membership. Other studies have used 
radial velocity or proper motion measurements to de- 


termine membership (e.g.. 

Dias et al. 2002 Nilakshi & 

Sagar|2002| IMermilliod et a 

1.120081 IFrinchaboy &: IVlajew- 

ski|2008|) , but these surveys have very few measurements 


and typically only cover the blue end of the MS or the 
cluster’s red giants. We therefore used the gi photome¬ 
try of HA08 for «16,800 objects with —0.5 < {g — i) < 


4 and combined it with the stars’ distances from M37’s 
center to determine membership probabilities (Pmem)- 
The full description of our Pmem calculations is in Ap¬ 
pendix!^ Stars with 0.2 < {g — i) < 3.6 were assigned a 
probability of being a single star (Pg), a binary star (Pb), 
or a field (non-member) star (Pf), so that Pg-I-Pb-I-Pf = 1 
for every star. We considered stars with Pmem = Ps 

See http://www.sdss.org/drl2/algorithms/ 

sdssUBVRITrcmsform 


+ Pb > 0.2 to be candidate cluster members. Mem¬ 
bers with Pb > Ps were labeled as likely binaries. The 
low Pmem threshold we used allows for a non-negligible 
number of contaminants into our cluster sample. How¬ 
ever, since we expect most field stars to be much older 
and, thus, inactive or significantly fainter in X rays, the 
contamination in our rotation-activity analysis of clus¬ 
ter stars will be minimal. In Figure we show the 
(i, g — i) CMD and we highlight single members, likely 
binary members, and field stars. 

We inspected objects with 0.1 < Pmem <0.2 that also 
have an X-ra y co unte rpar t with a hardness ratio <—0.3 
(see Sections 4.3|and [T4| ) as these could potentially be 
stars that are younger than the average contaminating 
field star, and hence potential cluster members. We 
found 12 objects that met these criteria, and we assigned 
them a Pmem value of 999. We use this Pmem value as a 
flag to highlight those stars that we consider to be cluster 
members based on additional criteria. 

For stars bluer than {g — i) = 0.2, we defined a separate 
CMD locus and calculated Pmem for those based only on 
distance from the cluster center and without distinguish¬ 
ing between Pg and Pb . This is because the MS is almost 
vertical in {g — i) for these stars, so that our standard 
approach cannot be applied. Objects considered clus¬ 
ter members based on these criteria are labeled as blue 
members. We show these as blue triangles in Figure]^ 

We did not apply our Pmem calculation to stars with¬ 
out HA08 gri photometry, namely, the 16 ME08 stars 
without an HA08 counterpart. Instead, we determined 
their membership by visually inspecting their location 
on the {i, g — i) CMD using t heir BV photometry trans¬ 
formed into gri (see Section 3.2 ) We also considered 
whether they matched to a soft (hardness ratio <—0.3) 
Chandra source. Six of these stars lie near or inside the 
locus of MS stars on the CMD and/or match a soft Chan¬ 
dra source. We assigned to these six either Pg = 999 or 
Pb = 999, based on their binary flag in ME08, and we 
indicate them with orange-filled squares in Figure 

We examined the literature for cluster me mbership in¬ 
forma tion for stars beyond the MS turn-off. jPias et al.| 
(20021 identified 40 cluster members near the turn-off 
point based on pro p er mo tion measurements, but |Frinch-| 


aboy & Majewski (20081 used radial velocity measure- 
ments to show that only three are M37 members. We 
found an HA08 match for one of these stars in our op¬ 
tical catalog. For the other two, we converted the avail¬ 
able BV photometry to gri using the transformations 
described in Section [3^ assigned them Pmem = 999, and 
added them to our optical catalog. 

We searched for surveys identifying white dwarfs and 
red giants in the field of M37. Usi ng a 5" tolerance ra- 

(p^l 

(SDOHl 


dius, four cluster white dwar fs in Kalirai et al. 
and 31 cluster red giants in |lVlerrnilliod et al. 
matched objects in our optical catalog; we assigned all 
these matched objects either Pg = 999 or Pb = 999, based 
on their binary flag in the literature. 


liod et al. 




in 

Mermil- 

Ni 

akshi & 


red giant and two massive stars matched X-ray so urce s 
that had no HA08 or ME08 counterpart (see Section 4.31. 
For these three cluster stars, we converted their PU pho¬ 
tometry to gri using the transformations described in 
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TABLE 5 

Catalog of Optical Objects: Basic Properties 


Position 


Properties 



Membership 

Derived Quantities 

Matched Sources 

a (J2000) 

& (J2000) 

9 

r 

i 

Prot 

Flag 

Ps 

Pr 

Mass 

T 

-^bol 

XID 

XID offset 

(°) 

(°) 

(mag) 

(mag) 

(mag) 

(days) 




{Mq) 

(days) 

(erg s-l) 


(") 

(2) 

(3) 

(4) 

(6) 

(8) 

(10) 

(11) 

(12) 

(13) 

(15) 

(16) 

(17) 

(19) 

(20) 

87.824771 

32.585669 

19.03 

17.96 

17.33 



0.36 

0 

0.78 

25.90 

7.80E-I-32 



87.887083 

32.658222 

17.48 

16.60 

16.12 

5.39 

2 

0 

0.04 

1.01 

17.11 

2.34E-I-33 

1 

1.2 

87.958613 

32.586522 

16.13 

15.51 

15.27 

3.92 

1 

0.51 

0.04 

1.19 

12.96 

6.37E-I-33 

40 

1.1 

87.966813 

32.596517 

22.42 

20.66 

19.43 

1.24 

1 

0.40 

0 

0.50 

48.35 

1.82E-I-32 

51 

0.6 

88.002083 

32.613389 

15.08 

14.15 

13.72 

1.89 

2 

0 

0 

1.57 


2.33E-I-34 

108 

0.1 


Note. — Only optical objects with either Proti Pmem > 0, or an X-ray counterpart are included in this list. This table is available in its 
entirety in the electronic edition of the ApJ. Five rows with some data columns are shown here for guidance regarding its form and content. 
Table [^describes all the columns in this table. 


Section 13.2 


_ assigned them either Pg = 999 or Pb = 999, 

based on their binary flag in the literature, and added 
them to our optical catalog. 

Our flnal membership catalog contains 1699 candidate 
M37 members, of which 265 are likely binaries. Twenty 
likely binary and 406 single members have Prot measure¬ 
ments. Table includes all optical objects to which we 
assigned a non-zero Pmem value. 

We note that there is a segregation of binary candi¬ 
dates toward the blue and red ends of the MS. That 
most stars at the red end of the MS are likely binaries 
may be a sign of mass segregation and tidal evapora¬ 
tion of the low-mass single stars, or it could be a sign 
of increased photometric uncertainty. Such mass segre¬ 
gation and preferential loss of low-mass single members 
has b een observed in man y other clusters, in cluding Prae- 
sepe (Pinfleld et al.|2003) and the Pleiades (Lodieu et al. 


20071, where the binary fraction increases in the lowest- 


mass bins in a way not seen for cool held stars (e.g., 
Duchene & Kraus 20131. Furthermore, the lack of bi- 
naries at 1.4 < [g — i) < 3.2 suggests that the binary 
sequence is not sufficiently well detected to be fit as a 
distinct component by our algorithm; the small number 
of over-luminous binary systems is likely being incorpo¬ 
rated into the fit for the held population, indicating one 
possible source of incompleteness. 

4. ANALYSIS 

4.1. Stellar Properties Derived from Optical Photometry 

4.1.1. Stellar Masses and the 500-Myr Mass-Period 
Distribution 

To estimate masses for all the stars in our optical cat¬ 
alog, we used the mass-absolute r magnitude (Mr) rela¬ 
tion of Kraus & Hillenbrand (20071. These authors com¬ 
piled a set of empirical observations from the literature 
to generate a complete list of derived masses and spectral 
energy distributions (SEDs) for B8-L0 stars, calibrating 
spectral types with stars in the 600-Myr-old Praesepe 


cluster (Kraus & Hillenbrand 


To obtain Mr we estimated the total absorption in r 


20071. 


(Ar) u sing the extinction tables by Schlafly & Finkbeiner 
(2011 [ ) assuming Ry = 3.1 and adopting a reddening of 


P(P - V) = 0.227±0.038 (HA08). We then adopted 
HA08’s distance to M37 of 1490±120 pc. Finally, we 
linearl y interpolated between the [Kraus fc Hillenbrand] 
(20071 Mr values to obtain masses. 

Our optical catalog includes stars rang ing from 0.16 
to 3.8 Mq. The top histogram in Figure 10 shows the 



-0.5 0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 

a - i 


Fig. 9. — CMD for all optical objects in the M37 field. Objects 
with Pmem > 0.2 are likely cluster members (blue, black, red, and 
orange symbols); objects with Pmem < 0.2 are likely field stars 
(gray points). Members with Pj, > Pg are flagged as likely binary 
members (red crosses). Blue members (i.e., those beyond the MS 
turn-off) have no single/binary flag (blue triangles). Orange-filled 
squares indicate ME08 stars without HA08 counterpart c lassified 
as member s. The faint blue members are white dwarfs l|Kalirai| 
jet al.|2005| |. The black arrow indicates the extinction vector. 

derived masses for cluster members, and the bottom his¬ 
togram includes only those members with X-ray coun¬ 
terparts in our Chandra catalog. We And X-ray emitting 
cluster members covering the entire surveyed mass range. 
We derived mass uncertainties from photometric and dis¬ 
tance uncertainties; we indicate typical uncertainty val¬ 
ues with gray diam onds and error bars in the bottom 
panel of Figure]^ the mass-period distribution for clus¬ 
ter stars. We distinguish between members with HA09 
clean or ME08 Prot and those with HA09 non-clean Prot 
in the upper panel (see Section [3T| ) , and between single 
and binary members with HA09 clean or ME08 Prot in 
the lower panel. 

There is a well-defined sequence of slow rotators in 
Figure [TT] that runs from M « 1.2 M©, Prot « 3 d, to 
M « OT^Mq, Prot ~ 20 d. Below this sequence is an¬ 
other distinguishable sequence composed mainly of rota¬ 
tors with non-clean Prot that are roughly half the Prot 
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Fig. 10.— Top: Stacked histogram of stellar masses for M37 
members. Likely binary members are indicated in red and blue 
members (which have no single/binary flag), in hashed blue. Bot¬ 
tom: Histogram of stellar masses for M37 members with an X-ray 
counterpart in our Chandra catalog. We do not derive the masses 
of the red giants and white dwarfs in our optical catalog and these 
are excluded from these histograms. 



of the slow rotators. This suggests that the second se¬ 
quence is composed of stars whose reported Prot are a 
harmonic of their true Prot • 

Although more than half of the cluster members with 
non-clean Prot he in or near the main mass-period se¬ 
quence and therefore display the same behavior as stars 
with clean Prot, we conclude that non-clean Prot add 
unnecessary ambiguity to our results. We therefore in¬ 
cluded non-clean rotators in our analysis only if they also 
had a ME08 Prot and the two Prot measurements did not 
differ by more than 20%, which is the typical vertical 
width of the main mass-period sequence. 

The mass-period distribution for M37 indicates that 
the more massive cluster members have spun down to 
the slow-rotating sequence, where Prot has a weak depen¬ 
dence on stellar mass. This breaks down around wO.T 
Mq; lower-mass stars span a range of 1 < Prot 20 
d. This is consistent with expectations based on M37’s 
slightly younger age tha n the Hyades and Pra esepe. In 
the 625-Myr-old Hyades (Perryman et al.|1998) and 600- 
Myr- old Praesepe clusters, this transition occu rs at «0.6 
Mq ( Agiieros et al.|[20TI Douglas et al.||2014 l. 

Cluster members with X-ray cou nter parts (filled sym¬ 
bols in the lower panel of Figure [TT| ) cover the entire 
mass-period plane, suggesting that A-ray detection is not 
just limited to the fast rotators (Prot ^ 1 d) and, hence, 
to the more active stars. Nine of the 18 fast rotators in 
our cluster sample have no X-ray counterpart. Eight of 
those nine rotators are at an off-axis angle 0 > 5' in the 
Chandra field-of-view, beyond which the Chandra PSF 
quickly degrades. Furthermore, based on their derived 
masses, four of them are likely to be early F stars, and 
the measured periods are probably due to pulsations and 
not to rotation. 



Fig. 11.— Top: Prot versus mass for M37 members. Our consoli¬ 
dated optical catalog has 426 cluster members with Prot: 296 have 
either HA09 clean Prot or ME08 Prot (blue stars), and 130 have 
HA09 non-clean Prot (gray crosses). A well-defined, slow-rotating 
sequence runs from M ^ 1.2 Mq, Prot ~ 3 d, to M ^ 0.4 Mq, 
Prot ~ 20 d. This sequence breaks down around a mass of ^0.7 
Mq-, lower-mass stars include both fast and slow rotators. Bot¬ 
tom: Same as the top panel, but for the 296 members with HA09 
clean Prot or ME08 Prot. Single members (283) are indicated with 
black circles, and likely binary members (l3) with red triangles. 
Stars with X-ray counterparts are shown as filled symbols (69 sin¬ 
gle members and seven likely binaries). Typical mass uncertainties, 
derived from photometric and distance uncertainties, are indicated 
along the top with error bars around gray diamonds. 

4.1.2. Bolometric Luminosities 


To obtain Lboi for the stars in our catalog, we once 
again adopted HA08’s distance to M37 of 1490T1 20 pc . 
We used the SED table of Kraus fc Hillenbrand| ( 20071 
to convert the Mr values obtained in Section |4. 1.1 1 into 
corresponding effective temperatures Peff- We then used 
Peff to find bolometric correction s BCr using the tables 
collected by Girardi et al.| (20041, which are tailored to 
the SDSS filter system. Using Mr and BCr, we then 
calculated bolometric magnitudes and luminosities. 


4.1.3. Rossby Numbers 

Stellar activity increases with increasing rotation rates 
only up to a saturation velocity, beyond which activity 
remains constant. This saturation point has been ob- 
served to have a dependence on stellar mass (see e.g., 
Pizzolato et al.|2003 1. Noyes et al. (19841 found that the 


Rossby number Kg = Brot/T, where r is the convective 
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turnover time, removes this mass dependence, thus re¬ 
vealing the general relationship between stellar rotation 
and activity. Ro serves as an indicator of the strength 
of a rotationally driven dynamo, which can be compared 
more evenly between stars of different masses. 

To estimate t for M37 me mbers, we used the e mpir- 
ical stellar mass-r relation of Wright et al. (2011 here¬ 
after Wll), who collected over 800 stars spanning a mass 
range of 0.09-1.36 Mq with both P^ot and Lx measure¬ 
ments. These authors fitted the Px/Tboi-7rot relation 
following a formulation equivalent to Equation and 
adopting the canonical slope f3 = —2.0. Wll reported 
r values for ten mass bins, and we used these to fit a 
second-order polynomial in the mass-r plane. The re¬ 
sulting mass-r relation is 

logr = 1.24-1.63 log(M/M 0 ) - 0.56 log^(M/M 0 ), 

with a correlation coefficient r = 0.989, and we used this 
relation to calculate r for cluster members within the 
0.09-1.36 Mq mass range. 72% of cluster members and 
99% of the cluster rotators fall within this mass range. 

Wll also reported r values for the same mass bins 
adopting /3 = —2.7, which is the slope they found in their 
analysis of a subset of 36 solar-type stars. We used the 
second-order polynomial that Wll fit to these r values 
in the mass-r plane to calculate r for our cluster stars, 
and found that these r values were 30±5% smaller than 
those based on /3 = —2.0. 

Using both sets of r values, we calculated Ro and ex¬ 
amined the Lx/L\^oi-Ro plane for our cluster sample, per¬ 
forming the parametrization analysis described in Sec¬ 
tion]^ twice. We found that the two (3 values we derived 
for our sample differed <3% and were well within one 
standard deviation (cr) of each other. The two satura¬ 
tion points Po.sat differed by 20 %, but were still within 
1(7 of each other. We therefore opted to use r values 
based on (3 = —2.0 for the rest of our analysis. 

4.2. Stellar Properties Derived from X rays 

The stellar magnetic field manifests itself above the 
photosphere and heats the plasm a in the stellar corona 
above 10® K ( Vaiana et al.|[l9^ ). X rays from MS stars 
originate in this high-temperature plasma. In such stars, 
the ratio of X-ray luminosity to bolometric luminosity 
Lx/Lhoi appears to decrease by several orders of magni- 
tude between zero-ag e MS stars and solar -age stars (e.g., 
Feigelson et al.||2004 iWright et al.||2010|. The Px/Aboi 
ratio thus allows us to compare tne contribution of X 
rays in a mass-independent way. 

To estimate the unabsorbed fluxes of our Chan¬ 
dra sources, we converted our soft-band fluxes using 
WebPIMMS^ and an APEC model with solar abun¬ 
dance and kT = 0.768 (log T = 6.95). Given that 95% 
of the X-ray sources with M37 counterparts have similar, 
soft X-ray spectra (wi th h ardness ratios between —0.4 
and —1.0; see Section 4.4), assuming the same plasma 


temperature for all is a reasonable approximation. We 
set the Galactic atomic neutral hydrogen column den¬ 
sity to 1.26x10^^ cm“^, derived using E{B — V) — 0.227 
from HA08, and obtained the fluxes for the ROSAT band 
(0.1-2.4 keV) in order to compare our results to those of 

http://cxc.harvard.edu/toolkit/pimms.jsp 
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Fig. 12.— Fraction of optical objects with X-ray counterparts 
per 10%-Pmem bin. 

previous studies. We then calculated Lx using a distance 
of 1490±120 pc (HA08). 

4.3. Matching the Optical and X-ray Catalogs 

We used the R:!80% confide nce X -ray positional uncer¬ 
tainties calculated in Section 12.21 to match our Chandra 
sources to optical objects from HA08 by using a tolerance 
radius 1.6 times the positional uncertainty of each X-ray 
source. We found that this factor provided a large num¬ 
ber of matches without introducing too many potentially 
false matches. We found 540 X-ray sources with HA08 
optical counterparts. Multiple matches occurred <15% 
of the time, and we selected the closest match as the true 
one in such cases. Fourteen more near-matches were in¬ 
spected individually and considered to be true matches. 

To match our X-ray sources to ME08 o bject s, we used 
the same method as described in Section |3.2[ where we 
described how a distortion in the absolute astrometry of 
ME08 objects led us to apply a liberal 10" tolerance ra¬ 
dius and then to inspect by eye the most likely true match 
using available gri and BV photometry as guidance. We 
found 37 X-ray sources with ME08 optical counterparts. 

For the remaining un matched X-ray sources, we found 
one red giant match in M ermilliod et al. ( 20 0 8t and two 

j 2002 |) (see 
.0 times the 
Neither of 


massive sta r matches in Nilakshi &: Sagar 


Section 3.3) also using a tolerance radius 


positional uncertainty of each X-ray source, 
these three optical objects has a HA08 or ME08 coun¬ 
terpart (two of them have V photometry outside the V 
range studied by ME08). 

Gombining all X-ray-to-optical matches, we have 561 
X-ray sources with optical counterparts, 278 of which 
have Pmem > 0.2: 195 single stars (69 with Prot) and 83 
binaries (7 with Prot)- Figure 12 shows the fraction of 


optical objects with X-ray counterparts as a function of 
Pmem- It is close to zero for field stars and then increases 
linearly to some intrinsic fraction of M37 stars. 

Two hundred and thirteen X-ray sources have no opti¬ 
cal counterparts. The log of the median spurious detec¬ 
tion probabilities P for matched and unmatched X-ray 
sources are —9.5 and —10.6, respectively. The median 
net counts in the full band for matched and unmatched 
X-ray sources are 17.2 and 25.9, respectively. Lastly, the 
median photon energy for matched and unmatched X-ray 
sources are 1.2 and 1.95 keV, respectively. 

4.4. X-ray Hardness Ratios and Spectral Information 
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Fig. 13.— Normalized stacked histogram of hardness ratios for 
X-ray sources matched to M37 members (white), X-ray sources 
matched to non-member objects (red), and X-ray sources with no 
optical counterpart (hashed blue). 


We obtained spectral information for our X-ray sources 
by computing a hardness ratio {HR). The HR is (hard 
counts—soft counts)/(hard counts-l-soft counts); we fol¬ 
lowed the definit ions of the hard and soft bands described 
in Section and used net counts. Soft spectra are 
characterized by HR < 0, and hard spectra, by HR > 0. 

An a p osteriori analysis of the matching done in Sec¬ 
tion |4.3| revealed that 96% of the 278 X-ray sources 
matched to optical cluster members have HR < —0.3, 
a strong indication of their stellar nature. Only two such 
sources have HR > 0, and both are fiagged as likely 
binaries. By contrast, 72% of the 283 X-ray sources 
matched to optical non-members (i.e., field stars) have 
HR < —0.3, and 12% have HR > 0. Furthermore, only 
32% of the 213 X-ray sources without any optical coun¬ 
terpart have HR < —0.3, and 45% have HR > 0. Fig¬ 
ure shows a histogram of HR values for X-ray sources 
matched to M37 members. X-ray sources matched to 
non-members, and X-ray sources with no optical coun¬ 
terpart. 

Only three X-ray sources matched to an optical clus¬ 
ter member have >150 net counts in the full band; the 
highest has 260. Thus, there are no cluster X-ray sources 
with enough counts for a meaningful analysis of their X- 
ray spectra or light curves. Four X-ray sources in our 
Chandra catalog have >1000 net counts. Two of these 
are matched to an optical non-member; the other two 
lack optical counterparts. We discuss these four sources 
in Appendix [B| 


5. RESULTS AND DISCUSSION 


Our analysis of Tx/.bboi versus Ro only included M37 
members (Pmem >0.2) that had either HA09 clean P^ot 
or that had both HA09 non-clean Pj-ot and ME08 Prot 
and the two Prot measurements differed by less than 20%. 
Figure 14 shows the resulting 51 single members (blue 
stars) and 4 likely binary members (red circles) that are 
included in the following fit. 

The conventional way to parametrize the Lx/Lboi~Ro 
relations hip is with a fiat region connected to a power- 


law (e .g., Randich 2000 Wright et al. 2011| Douglas et al. 


20141: 


\ CPf if Po > Po,,at 


( 2 ) 


where (Tx/Aboi)sat activity saturation level, Ro,sat 

is the turnover point, /3 is the power-law index character¬ 
izing the unsaturated regime, and C is a constant. We 
fit this model to the 55 members described above. 

We follo wed a similar fi t ting p rocedure to that de¬ 
scribed in Douglas et al. (20141. We converted the 
model in Equation into log-space for the fit, and as¬ 
sumed fiat priors over all three parameters. We used the 
open-source Markov-chain Monte Carlo (M CMCl pack¬ 
age emceq^ (Foreman-Mackey et al. 20131 to carry out 
the fit. It Duilds posterior probab ility distributions for 
each parameter (shown in Figure 151 by performing a 
random walk in parameter space. 

Two hundred models drawn at random from the pos- 
teri or p robability distributions are shown in gray in Fig¬ 
ure |14| The most probable model, i.e., the maximum 
a posteriori model, is indicated by a solid black line. 
The resulting best-fit parameters are log(Px/Aboi)sat = 
-3.06 ± 0.04, Po,sat = 0.09 ± 0.01, and /3 = -2.03ti];J^. 


The parameter values correspond to the 50**' quantile, 
and the uncertainties correspond to the 16*** and the 84**^ 
quantiles (for consistency with Icr Gaussian uncertain¬ 
ties). 


5.1. The Saturated Regime 

Several hypotheses aim to explain the decoupling of ro¬ 
tation and activity seen in the saturated regime. These 
include d ynamo saturation (Gilman||1983 Vilhu & Wal- 
ter[1987 ), saturation of the filling factor of active regions 
on the stellar surface ( Vilhu|19 84p, and cent rifugal strip¬ 
ping of the corona (Jardine & Unruh 19991. However, a 
satisfactory explanation remains elusive (see e.g., Wll). 

We detect this saturation phenomenon for stars with 
Ro < 0.09 ± 0.01. This turnover point is simi¬ 
lar to Ro,sat = 0.13 ± 0.02 found by Wll for their 
Einstein/ROSAT/XMM heterogeneous sample of 824 
field and cluster stars collected from the litera ture, and 
i t is sl ightly smaller than Ro,sat ~ 0.16 found by|Randich 
(20001 for her heterogeneous ROS AT sample of held and 
cluster stars. Our saturation level log(Lx/iboi)sat = 
—3.06 ± 0.04 agrees with log(Lx/iboi)sat = —3.13 ± 0.22 
found by W1 1. It is slightly lo wer than log(Lx/Aboi)sat ~ 
—3 found by Randich ( |2000[ |. 


5.2. The Unsaturated Regime 

The increase in stellar activity with increasing rota¬ 
tion (smaller Ro) is believed to i ndicate the heating of 
the corona by th e stellar dynamo (Pallavicini et al.|1981 


Pizzolato et al.|[2003 1. One of the challenges ot the past 


tew decade s has beeii to characteriz e empirically this re¬ 
lationship. Pallavicini et al. (1981) first discovered this 
behavior in an Einstein saniiple of 17 late-type stars of 
luminosity classes II to V, and fou nd a power-law in- 
dex /3 = —1.9 ± 0.5. More recently, iGiidel et al.| (|1997 1 
found (3 = —2.64 ± 0.12 for a ROSA'T/ASCA sample of 
12 solar-type stars of ages 0.07 to 9 Gyr, and | Randich] 
([2000l found /3 = -2.10 ± 0.09 for her ROSAT sam¬ 
ple. Finally, Wll found /3 = —2.18 ± 0.16 for their het¬ 
erogeneous sample of cluster and field stars. Wll also 
analyzed a subset of 36 solar-type stars detected only 
with ROSAT to eliminate X-ray luminosity biases. They 


http: //dcLQ. iel. fm/emcee/current/ 
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Fig. 14. — X-ray to bolometric luminosity ratio as a function of 
Rossby number Ro for cluster members. Likely binary members 
are indicated with red circles. The solid black line is the maximum 
a posteriori fit from the MCMC algorithm and the gray lines are 
200 random samples from the posterior probability distributions . 
The dashed line indicates the fit found bylPa llavicini et al. 1 1981 1, 
and the dotted line indicates the fit found by |wrignt et al. (^UTTl. 
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Fig. 15.— Marginalized posterior probability distributions from 
the MCMC analysis using emcee. The parameter yalues of the 
a posteriori model are the peaks of the one-dimensional distribu¬ 
tions; the yertical dashed lines approximate the median and 68- 
percentiles. The two-dimensional distributions illustrate coyari- 
ances between parameters; the contour lines approximate the 68- 
and 95-percentile of the distributions. 


found /3 = —2.70±0.13 for this sample. We include the 
power-law relation for th e unsaturated regime of b oth the 
latter sample and of the Pallavicini et al. (19811 sample 
in Figure [HI 

For our sample of unsaturated, co-eval M37 stars, 
which cover a mass range of 0.55-1.28 Mf^, the index 


/3 = —2.03tli4 agrees with that by 


Pallavicini et al. 


(19811 j3 = —1.9 ± 0.5. It differs from the Wll result for 


the unbiased sample of solar-type stars by 3 ct. That /3 
found for our sample of ~500 Myr stars does not agree 
with some of the /3 values in the literature suggests that 
previous studies that include stars with a range of ages 
may be missing interesting age-dependent effects in the 
behavior of stellar activity. It may also be evidence that 
the unsaturated regime is better described by a broken 


powe r law instead of a single power law (e.g.. Petit et al. 
[Donati fc Landstreet||200^ [Saar fc Brandenburg 


2008 


19991 Unfortunately, our sample of M37 does not ex¬ 
tend to large enough Ro values to fully test this idea. 

Did our Chandra observation or the P^ot surveys of 
ME08 and HA09 bias our sample toward those stars with 
higher Lx and/or faster Prot? We may not have identi¬ 
fied the faintest X-ray sources in the cluster because of 
their very low counts and/or large 9, making their de¬ 
tection much harder. Likewise, ME08 and HA09 may 
not have measured Prot for the faintest cluster members, 
thus excluding members with possibly the longest Prot- 

Our X-ray catalog contains 18 sources, 15 of which 
are cluster members, with log(Lx/Lboi) < —6.0, while 
our optical catalog contains only six objects, three of 
which are members, with Ro > 1.0. Indeed, the light 
curves used by both ME08 and HA09 span only 30 days, 
suggesting that lack of slow rotators may be the limiting 
factor in determining the behavior of unsaturated stars 
at the high-Po end. 

Still, our sample of «500 Myr M37 members with both 
Prot and X-ray measurements spans a wide set of stellar 
properties, as it covers masses ranging from 0.4 to 1.3 
Mq, Prot ranging from 0.4 to 12.8 d, and Lx ranging 
from 10^® to 10®° ® erg s“^. Our results make M37 a 
new benchmark open cluster for calibrating the rotation- 
activity relation at ages of «500 Myr. 


5.3. The Dependence of X-ray Activity on Age 

We compared Lx for 64 M37 solar analogs (masses 
0.8-1.2 Mq) with X-ray detections of 62 stars in Alpha 
Persei (a-Per, Randich et al. 1996 Prosser et al. 1996 


Prosser fe Randichl 119981 [Prosser et al.||19981 |Pil Itten 


et aL|2013p, 62 Plemds ([Micela et al.|1990| 


Stauffer et al. 

2003|, 28 

19941119951), and five stars in NGC 


TUM 


iVliceia et al. 


Hya ds (Stern et al. 


199' 


752 (Giardino et al. 


MThiriggs fc 


2008 ) within the same mass range. 


The a dopted ages are 85±10 Myr for g-Per ([Lodieu et ah 


20121), 130T20 M yr for the Pleiades ( RarrLcT^^^ras^ 
cues et ai.||2004), 48 5T28 Myr for M3 y (H A08), 625±50 
IVlyr for the HyLes (Perryman et al. ]I9^ , and 1.6±0.2 
Gyr for NGG 752 (Agiieros et al. iii prep.). Figure 16 
shows the data for these five clusters, where the boxes 
extend from the lower to the upper quartile, the whiskers 
extend to cover the entire range, the red line inside the 
boxes indicates the median, and the blue stars indicate 
the mean. The m aximum and minimum Lx of the Sun 
(Peres et al. 2000) is indicated by black whiskers at 4.5 
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Fig. 16.— X-ray luminosity of solar-like members (mass range of 
0.8-1.2 Mq) of five open clusters as a function of age. The adopted 
ages are 85dil0 Myr for a-Per, 130di20 Myr for the Pleiades, 
485±28 Myr for M37, 625±50 Myr for the Hyades, and 1.6±0.2 
Gyr for NGC 752. The data includes 62 a-Per stars, 62 Pleiades 
stars, 64 M37 stars, 28 Hyades stars, and 5 NGC 752 stars. For 
each cluster, the box extends from the lower to the upper quartile, 
and the whiskers cover the entire Lx range. The median Lx for 
each cluster is indicated by a red line inside the box, and the mean 
by a blue star. We in dicate the minimum and maximum of the 
Sun (|Peres et al.||2000|l with black whiskers at 4.5 Gyr. A linear 
regression analysis m log space (excluding the Sun) suggests a re¬ 
lationship between X-ray luminosity and age of Lx oc 
and we indicate it with a dashed line. 

Gyr. 

We performed a linear regression analysis in log space 
using the median Lx values of the five clusters to deter¬ 
mine the dependence of Lx on age. We found that Lx oc 
t“i-23±o.i6^ with a correlation coefficient r = —0.98. This 
result, which corresponds to the black dashed line in Fig¬ 
ure is close to the expected power-law exponent value 
of —1 for stars with ages ,>100 Myr. It is within la of 


the power-law exponent — 1.5 ^q2 found by Maggio et al. 


(19871 using a set of Pleiades, Hyades, and Orion Ic solar¬ 
like stars and their observatio ns o f nearby solar analogs . 
Similarly, Giidel et al. (19971 and Giardino et al. (20081 
found an exponent of —1.5. Plugging our (3 = — 2.037ou4 
into the relation Lx oc we find that a must lie in the 
range 0.49-0.75, with a = 0.60 as the most likely value. 
This agrees very well w ith the relation Vm t oc ® for 


solar- l ike sta rs found by Skumanich (19721 and Meibom 

eFaTlp^. 


6 . SUMMARY 

We have assembled an extensive sample of stars with 
X-ray and P^t measurements in the M37 open cluster to 
study the stellar rotation-activity relation at «500 Myr. 
We observed the center of the cluster with Chandra and 
detected «770 high-confidence X-ray sources. We also 
consolidate d two optical surveys of M37 with Prot mea - 


surements (Messina et al. 2008 
obtaining a merged sample of 


Hartman et al. 20091, 
650 rotators. We calcu- 


lated new membership probabilities for optically detected 
objects in the M37 field, and complemented the resulting 
catalog of cluster members with massive stars, red giants, 
and white dwarfs identified as members in the literature. 
By matching our X-ray and optical catalogs, we gener¬ 
ated a sample of 76 solar- and late-type members of M37 
with X-ray and Prot measurements. 

We used soft-band X-ray fluxes extracted from our 
Chandra observations and the cluster distance of f«1.5 
kpc to derive Lx for all X-ray sources. We then used gri 
photom etry and reddening values from [Hartman et al.| 
(2008a) to derive stellar masses and Lboi for all clus- 
ter members. Finally, we calculated r values using the 


Wright et al. (2011) empirical model as a function of 
mass, thereby obtaining Ro for our stars. 

We then examined the «500 Myr mass-rotation and 
rotation-activity relations. M37’s mass-period distri¬ 
bution indicates that the cluster’s more massive stars 
have spun down to a single-valued relation. There is 
a well-defined sequence of slow rotators that runs from 
M w 1.2 Mq, Prot w 3 d, to M w 0.4 M©, P^t ~ 20 d. 
This single-valued relation breaks down around a mass 
of 0.7 M©; lower-mass stars span a range from fast (Prot 
~ 1 d) to slow («20 d) rotators. This is consistent with 
expectations based on M37’s slightly younger age than 
the Hyades and Praesepe. In these ftiOO O-Myr-old clus¬ 
ters, this transition occur s at «0.6 M© (Agiieros et al. 
2011 Douglas et al.||2014 l. 


TiT examining the dependence of Lx/Lboi on Ro, we 
identified a saturated and an unsaturated regime. All 
rotators with Ro < 0.09 ± 0.01 (a turnover point consis¬ 
tent with that found in the literature for other rotation- 
activity samples) are saturated, and converge to an ac¬ 
tivity level of log(Lx/Lboi)= —3.06 ± 0.04, close to the 
canonical value of —3 and similar to recent values found 
for heterogeneous samples of field and cluster stars. 

In the unsaturated regime, faster rotators (smaller Ro) 
have increasing levels of coronal activity (larger Lx/Lboi) 
following a power law of index /3 = —2. 03tn'ili in agree¬ 
ment w ith the /I = —1.9 ± 0.5 found by [Pallavicini et ah 
(|1981t. Interestingl y, however, this j3 differs from the 
Wright et ai.| (|2011|) result for their unbiased sample of 
solar-type stars by 'Sa. This suggests that previous stud¬ 
ies that include stars with a range of ages may be missing 
interesting age-dependent effects in the behavior of stel¬ 
lar activity, or that a broken power law may describe 
better the unsaturated dependence of activity on rota¬ 
tion. 

A comparison of solar-type X-ray emitters in M37 and 
similar stars in the a-Per, Pleiades, Hyades, and NGG 
752 open clusters, covering an age range of 0.85-1.6 Gyr, 
shows that Lx oc i-i- 23 ±o.i 6 ^ suggesting that the sur¬ 
face rotational velocity in solar-type stars declin es as 


Vrnt oc t in a g reement with previous studies Sku¬ 


manich (e.g., 1972); Meibom et al. 


r revio 
(e.g.. 


2011 ). 


By probing the stellar rotation-activity relation with 
our homogeneous, co-eval sample of M37 members, we 
are helping to calibrate empirically the ARAR at «500 
Myr. More importantly, the wide coverage in stellar 
masses, rotation periods, and X-ray luminosities of our 
sample of cluster members puts M37 in a position to be¬ 
come a new benchmark middle-aged open cluster. 
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Fig. Al.— CMD for all optical objects in the M37 field. The dashed line indicates the MS defined in this paper by fitting a polynomial 
to a subset of stars close to the cluster center t hat made the MS visually obvious. The blue dotted line indicates the color-magnitude 
sequence derived by [Kraus &: Hillenbrand] ||2007[l for Praesepe. 


APPENDIX 

A. CALCULATION OF MEMBERSHIP PROBABILITIES 

ME08 and HA08 relied solely on distance from a fiducial isochrone in a CMD to determine membersh ip for M37. 
Other studies have used radial velocity or proper motion measurements to determ ine membership (e.g., Dias et al. 
2002 Nilakshi & Sagar||2002 Mermilliod et al. |2008 Frinchaboy & Majewski||2008 ), but these surveys have very tew 
measurements and typically only cover the blue end oi the MS or the cluster’s red giahts. Since kinematic measurements 
of objects in the field of M37 are scarce, we decided to use the spatial distribution of stars, in addition to their location 
in a CMD, to determine membership probabilities. 

A.l. Description of the Method 

First, we binned all the stars in the M37 field by {g — i) color, which allows us to keep both single and binary 
cluster stars within the same bins. We then defined a MS in the {i, g — i) CMD by fitting a polynomial to a subset 
of stars close to the cluster center that made the MS visually obvious. Figure [XT] shows the full set of optical objects 
in the M37 field, and we indicate with a dashed line the MS we defined. We also include the color-magnitude relation 
derived by Kraus & Hillenbrand (20071 for Praesepe members. The slight discrepancy between the two sequences is 
expected due to the difierences in filter systems and in metallicity, age, and reddening. Next, we built probability 
distribution functions (PDFs) for the single, binary, and field star populations. The two quantities that we used for 
fitting in each {g — i) bin were the radius r from the cluster center and height m = i — zms above or below the MS in 
the (i, g — i) CMD. By thinking in terms of height above/below the sequence (i.e., fitting and subtracting the MS), we 
worked in terms of a quan tity for which the distribution of contaminants is roughly constant across the bin at a given 

shows the distribution of m for stars in the M37 field in terms of {g — i) color and radius 


height/depth. Figure 


A2 


Stars are color-coded using the membership probabilities obtained after performing the fit 


from the cluster center, 
described next. 

The parameters that we fitted in the PDFs were: 

1. Three normalization constants: the number of single cluster stars (W)) of binary cluster stars (Ag), and of field 
stars (Nf). Clearly, As-|-Ab-l-Af=Atot- Since Atot is set, there were really only two parameters being fitted. 

2. One spatial distribution parameter: the exponential scale radius of the cluster tq. We assumed a fiat spatial 
distribution of field stars, which is defined by the size of the field and Af and, hence, does not introduce any 
additional fitting parameters. 


3. Two parameters setting the intrinsic width of the MS in terms of a of m: one around the single star sequence 
(us) and one around the binary star sequence (ub)- 
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TABLE A1 

Cluster Fit Parameters for M37 


SpType 

N{ 

Ns 

Nb 

ro 

(deg) 

a 

ms 

(mag) 

O-s 

mb 

(mag) 

o-b 

F5.0-G8.0 

137 

165 

23 

0.044 

-0.238 

0.010 

0.181 

0.636 

0.260 

G8.0-K5.0 

915 

162 

87 

0.043 

-0.395 

0.018 

0.103 

0.373 

0.306 

K5.0-M0.0 

446 

156 

0 

0.054 

-0.309 

0.027 

0.068 

0.357 

0.347 

M0.0-M2.5 

468 

153 

0 

0.053 

-0.348 

0.017 

0.084 

0.535 

0.448 

M2.5-M4.5 

647 

96 

0 

0.050 

-0.394 

0.050 

0.119 

0.455 

0.058 

M4.5-M6.0 

940 

18 

53 

0.046 

-0.455 

0.017 

0.064 

0.244 

0.192 


4. Two parameters setting the median height of the MS: one for the single star sequence (wg) and one for the binary 
star sequence (mb). 

5. One parameter a setting the (linear) slope of the field star brightness distribution within the bin, dN/d(Am). 
By defining the field star distribution in terms of height above/below the cluster sequence rather than in terms 
of i, we linearize what is otherwise a two-dimensional function of i and {g — i). 

The PDFs for the cluster single star population, cluster binary star population, and field population are therefore 
given, in units of stars per unit area per unit magnitude, by 


1 r 1 

g(m,r) = iVg X [-—2 exp(- )] x [ 

ZTrrn rn 




exp(- 


(m - rrisY 

2^r“ 


(Al) 


0 b(m,r) = iV 6 X [^;^exp(-^)] x exp(- )], and 


2(jI 


(A2) 


0 f (m, r) = [Atot - Ns - A’b] x [-^] x [a(m - 7 ) + ^] 

^tot 4 Z 


(A3) 


Note that the final term in Equation |A3[ which treats the magnitude dependence linearly, contains constants of 
normalization that change depending on the interval over which the field population is being fitted. This calculation 
assumes that the lower limit is —0.75 mag and the upper limit is -1-1.25 mag. 

A.2. Application of the Method to M37 

We considered all stars in a color bin that fell within 0.75 mag below or 1.25 mag above the MS. We fit for the values 
for all of these parameters in each color bin using a random walk around an initial set of guessed parameter values to 
find a minimum in space. We used 

0.2 <{g-i)< 0.7 (F5-G8), 

0.7 <{g-i)< 1.4 (G8-K5), 

lA< {g-i)<2.l (K5-M0), 

2.1 <{g-i)< 2.8 (M0-M2.5), 

2.8 < {g- i) < 3.2 (M2.5-M4.5), and 

3.2 <{g-i)< 3.6 (M4.5-M6) 

as the color bins covering the range of the discernible MS. We then calculated three probabilities to each star: one 
of being an isolated member (Pg), one of be ing likely a binary (Pb), and one of bei ng a fiel d sta r (P f). Note that 
Ps + Pb + Pf = 1.0 for every star. Table [AT] summarizes the parameters of Equations |A 1 [ |A 2 [ and |A3| obtained with 
the fit, and Figure [A3| shows a histogra m of the number of candidates as a function of Pmem for each color bin. 

We observe in the left panel of Figure |A2| a segregation of likely binaries toward the blue and red ends of the MS. 
The fact that most stars in the red end of the sequence are considered binaries illustrates a possible sign of mass 
segregation and tidal evaporation of the low-mass single stars, so that only the more massive binaries are left. This 
mass segr egation and preferen tial loss of low-m ass single members has been seen in many other clusters, including 
Praesepe (Pinfield et al.|2003l and the Pleiades (Lod ieu et al.||2007t where th e binary fraction increases in the lowest- 
mass bins in a way not seen tor cool field stars (e.g., Duchene & Kraus 2013). The lack of binaries at 1.4 < {g — i) < 

3.2 (K5-M4.5) suggests that the binary sequence is not sufficiently well detected to be fit as a distinct component by 
our algorithm. The small number of over-luminous binary systems is likely being incorporated into the fit for the field 
population, indicating one possible source of incompleteness. 

Lastly, we treated stars blue of the turn-off point (at (g-i) « 0.2) differently because the MS there is almost vertical 
in {g — i) and, thus, our standard approach cannot be applied. Instead of having height m above or below the MS 
as one of the membership criteria, we considered as a single group all stars within the following GMD locus (with 
no further consideration of GMD position or single/binary flag): —0.3 < {g — i) < 0.2 and brighter than 0.75 mag 
below the MS, and 0.2 < {g — i) < 0.5 and brighter than 1.25 mag above the MS. This locus catches stars blue of 
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g-i Radius from cluster center (degrees) 

Fig. A2.— Left — Distance m above/below the MS of M37 in the (f,(g — i)) CMD as a function of {g — i) color for stars in the M37 
field. We color-code stars using our calculated Pmem: field stars (Pmem < 0.2) in gray, single members (Pmem > 0.2 and Pa > Pj,) in 
black, and likely binary members (Pmem > 0.2 and P^, > Pa) in red. Vertical dashed lines indicate the limits of the color bins used in the 
analysis. Right — Same as left panel, with distance from the MS of M37 as a function of radius from cluster center. Binaries appear to 
be concentrated toward the center of the cluster. 
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Fig. A3.— Normalized histograms of the number of cluster candidates (Pmem > 0) as a function of Pmem for each color bin. The total 
number of stars N in each color bin is shown. The vertical dotted line indicates the value of Pmem = 0.2 that we use in our analysis of 
M37 to separate cluster member from non-members. 


the turn-off and those moving across the Hertzsprung gap. We then fitted using the spatial term as before to factor 
into the membership probability. Applying this approach, we found that of the 198 stars in the locus, 171 have 


> 0.2, 139 have Pmem > 0 -5, and 74 have P„ 


0.8. The scale radius is 0°.038, which is smaller than that of the 


redder bins (see Table Al I. This is consistent with the suggestion of mass segregation seen in the rest of the cluster 
members. 


A.3. Testing the Method with Praesepe 
To test the robustness of the membe rship probability method de scribed above, we implemented it on stars in 


and around the Praesepe open cluster. Kraus & Hillenbrand (2007) calculated Praesepe cluster Pmem using both 


photometric and proper motion data to identify wlOOO members from among the several million stars they examined. 

To implement our membership calculation, which excludes proper motion information, we limited our sample of 
Praesepe point sources to a circular locus in the proper motion diagram centered around the mean cluster proper 
motion. This was needed to obtain a sample set with a members-to-contaminants ratio similar to that of the M37 
field. We also limited our sample to objects in the range 2.8 < (g — i) < 3.2 (spectral types «M2.6-M4.3,) which gives 
us enough cluster members to compare. 
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mem 


Fig. A4.— Comparison between two different sets of Praesepe cluster Pmem cal culations: one using CMD distance, radial distance, and 
proper motion information calculated by Kraus Hillenbrand i 20071, and another one using only CMD and radial distances 

(PmenP'^^)^ calculated here. As expected, the overwhelming majority of stars have comparable Pmem: likely members are in the upper 
right locus, and non-members in the lower left locus. The subset of stars with high PmenP"*”^ (upper left) are likely part of 

the wings of the actual member distribution. 


The resu lts of our own CMD+radial-distance-based membership probability calculations are s hown in 

Figure A4 where we compare them against proper-motion-based calculations by Kraus fc Hillenbrand] (20071 
for the same stars. We found 296 objects with > 0.6 (232 single stars and 64 binaries); these are high- 

confidence Praesepe members. Kraus & Hillenbrand (20071 found 310 members in the same photometric range. 

The expectation that stars with high should also have high app ears to hold in practice. Likely 

members flagged as such by both calculations are in the upper right of Figure |A4| while likely field stars are in the 
lower left. Sources in the upper left (high low may be either cluster members on the wings of the 

photometric distribution or interlopers slightly in front or behind the cluster. The progressive way in which this 
subset of members smoothly meets the member locus in the upper right convinced us that these are most likely part 
of the wings of the actual member distribution. Objects in the lower right (low P^^, high are probably 

non-member dwarfs in the same spatial volume considered in the analysis. In summary, we deem our membership 
probability calculation to be robust. 

B. SPECTRAL AND TIMING ANALYSIS FOR SOME HIGH-COUNT X-RAY SOURCES 

We performed spectral and timing analysis on the four highest-count sources in our Chandra catalog: XID 65 (1267 
total source counts in the full band), XID 108 (3097), XID 144 (439 8), a nd XID 1004 (5814). Only XID 65 and XID 
108 have counterparts on our catalog of optical objects (see Section 3.2 ) but neith er one is a likely cl uster member. 


Source XID 108 is matched to the eclipsing binary of W UMa type KVll, which Kang et al. 


(2007) characterized 


as having a period of 1.89 days and a V magnitude amplitude of 0.15 mag. On the other hand, the only knowledge 
we have of the optical counterpart to source XID 65 is its gri photometry. Table |B1| shows the characteristics of the 
optical counterparts to the high-count X-ray sources, as well as the hardness ratio of the matched X-ray source. 

We used the QIAO tool Sherpc^^ to fit each spectra with three different xspec models using solar abundances 
by Wilms et al. (2000)): one-temperature (IT) and two-temperature (2T) APEC models, and a power law model. 
APEC models are used to model the X-ray emission of diffuse gas in collisional equilibrium, characteristic of stellar 
coronae. The power-law model is a simple photon power law of the form A{E) = KE~^, where AT is a normalization 
constant. W e use d the ISM absorption model tbabs using photoelectric cross-section from |Balucinska-Church fc| 


McCammon (1992) to account for extinction by atomic neutral hydrogen. Since the distance to these X-ray sources 
is undetermmed, we left the atomic neutral hydrogen column density (Vjj) as a free parameter. For each source, 
we fitted the spectra of the five individual observations simultaneously using the chi-squared (x^) statistic with the 
Gehrels variance function and Sherpa’s levmar optimization method. We obtained la or 68% confidence intervals 
of all free parameters by computing the covariance matrices. In Figure B1 we compar e th e binned spectra (15 to 25 
counts per bin) of the sources against a best fit IT-APEC model. We report in Table B2 all acceptable fits for each 
X-ray source. 


We only subtracted background counts on source XID 1004, as this source is located on the back-illuminated, more 
background sensitive ACIS-S3 chip. The ratio between background and net count rates for XID 1004 is «40%. For 
the other three sources, their location on the front-illuminated ACIS-I chips guarantees that background noise remains 
low. The ratio between background and net count rates is 8.5% for XID 65, 3.6% for XID 108, and 4.2% for XID 144. 

We used AE’s tools to extract adaptively smoothed light curves of the four sources. These tools also perform a 
test of time variability by comparing a uniform firrx model to the distribution of source event time stamps and then 

http://cxc.harvard.edu/sherpa/ 
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TABLE B1 

Properties of Optical Counterparts to High-Count X-ray Sources 


XID 

Opt. ID 

9 

(mag) 

r 

(mag) 

i 

(mag) 

V 

(mag) 

B-V 

(mag) 

Fx/FopC 

HR 

65 

240271 

21.53 

20.76 

20.46 

21.06'" 


2.88 

-0.12 

108 

220858 

15.08 

14.15 

13.72 

14.56 

0.66 

0.01 

-0.72 


^ Ratio deri ved using the absorb ed X-ray flux in the 0.5—7 keV band and the flux ratio 
deflnition by |Stocke et ar| ) |l99l1 l. 

^ Derived using the gr-to-V transformation by [jester et al^ | |2005[ . 


TABLE B2 

Spectral Fits for the Highest-Count X-ray Sources in our Chandra Catalog 


XID 

Model 

kTi 

(keV) 

kT2 

(keV) 

r 

Nh 

(cm-2) 

xlO'^i 

Flux 

(erg cm“^ s“'^) 
xlO”"'' 

xl/d.o.f. 

Note 

65 

IT 

10.57±2.11 



3.16±0.60 

5.44±0.27 

0.43/79 



PL 



1.56±0.12 

3.77±1.01 

5.59±0.32 

0.44/79 


108 

IT 

1.54±0.05 



0.45±0.45 

5.19±0.28 

1.50/140 

Flare subtracted 


PL 



3.59±0.17 

3.76±0.48 

13.31±1.72 

1.19/140 

Flare subtracted 


2T 

0.79±0.10 

3.77±1.91 


6.25±1.13 

41.00±5.69 

0.69/25 

Flare only 


PL 



2.60±0.23 

3.04±1.00 

30.90±4.28 

0.69/26 

Flare only 

144 

IT 

6.97±0.75 



4.47±0.33 

19.03±0.49 

1.32/156 



PL 



1.79±0.07 

5.83±0.56 

20.51±0.72 

1.31/156 


1004 

IT 

6.84±0.68 



4.89±0.35 

27.18±0.73 

0.91/344 

Background subtracted 


2T 

6.26±0.67 

0.34±0.18 


5.94±0.74 

29.29±1.12 

0.91/343 

Background subtracted 


PL 



1.78±0.07 

6.10±0.48 

29.92±0.90 

0.91/344 

Background subtracted 


Note. — Results of all acceptable spectral fits for our four highest-count X-ray sources. For the one-temperature (IT) and 
two-temperature (2T) APEC models we give the temperature components kT\ and fcT 2 . For the power law (PL) model we give the 
photon index F. All fits included the xspec model tbabs to account for ISM extinction, and we give the best fit atomic hydrogen 
column density Ah- We also show the derived unabsorbed flux in the 0.5—7 keV band. Finally, we give for each fit its reduced 
chi-square (xt) degrees of freedom (d.o.f.). 


computing a 1-sample Kolmogorov-Smirnov statistic Pks- A statistic value of Pks >0.05 means there is no evidence 
for variability, while Pks < 0.005 indicates that the source is definitely variable. Only XID65 had Pks > 0.05, 
suggesting that it is the only one of the four sources with no evidence of variability. In Figure [B2] we put in sequence 
the light curves of our five Chandra observations for each source, and we show Pks for each source. The vertical dotted 
lines indicate the cutoffs between observations, and the time between observations is not included. Note that source 
XID 1004 was not in the last observation (ObsID 13657). 

We used AE’s ae_pileup_screening tool to check for potential photon pile-up in these four high-count sources. 
The procedure calculates on the individual observations the number of counts per ACIS frame using 3x3 pixel islands 
within the PSF of the source. We found that the highest 3x3 rate was 0.02 counts per frame—corresponding to a ~1% 
pile-up fraction—for XID 108, and 0.01 counts per frame for XID 144. These pile-up fractions are below the suggested 
0.05 counts-per-frame value above which pile-up reconstruction is warranted. Sources XID 65 and 1004 had a 3x3 rate 
of zero. 

There is a multi-standard deviation flare in source XID 108, as evident on its light curve in Figure [B^ Including this 
flare in the spectra of the source led to unacceptable spectral fits {xt > 2) on all the models we tried. We thus fi tted 
the spectrum of the flare separately from the rest of the spectrum of XID 108, and we included both fits in Table |B^ 

All four spectra are well-fitted using either a IT-APEC model or a power-law model. Only source XID 1004 and the 
flare of XID 108 had a 2T-APEC model fit that resulted in two kT temperatures significantly above zero. Finally, the 
A^h values derived by the tbabs mod el fo r the four X-ray sources are larger than the 1.26x10^^ cm“^ for M37 (derived 
from optical reddening, see Section 4.21, suggesting that these four high-count sources may be background sources. 
The only exception is one of the fits for XID 108: A^h from the IT-APEC model fit suggests that it is a foreground 
object, whereas A^h from the power-law model fit suggests that it is a background object. The web-based FTOOL 


ibiec 

n/fij indicates a total Galactic A^h of «4xl0^^ cm“^ in the direction of M37 based on the Leiden/Argentine/Bonn 
survey of Galactic atomic hydrogen. This means that sources XID 144 and XID 1004, having A^h larger than that, are 
most likely extragalactic sources. 
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http://heasarc.gsf c.nasa.gov/cgi-bin/Tools/w3nh/wSnh. 
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Fig. B1. — Spectral fits for the brightest sources in our Chandra catalog. The IT-APEC model fit for each spectrum (solid gray line) is 
drawn over the binned (15-25 counts per bin) data points. The error bars of the binned data points are too small to show. The parameter 
values for the fit using the IT-APEC model are included in each panel: temperature of IT-APEC model (kT), atomic neutral hydrogen 
column density {Nn), reduced chi-squared statistic (xj), and degrees of freedom (d.o.f.). The residuals of the fit are shown at the bottom 
of each panel, normalized to the sou rce counts. The source spectrum and spectral fit of XID 108 exclude events during a multi-standard 
deviation flare, as seen in Figure [B^ 



Time (ks) Time (ks) 

Fig. B2.— Adaptively smoothed light curves of the brightest sources in our Chandra catalog in the 0.5-7 keV energy range. Background 
counts are not subtracted. Vertical dotted lines represent the cutoffs between our different Chandra observations. The light curves from 
the five observations are in chronological order from left to right, and the time between observations is ignored. Note that source XID 
1004 was not in the last observation (ObsID 13657). The Pks statistic of variability is shown for each source. Pks > 0.05 indicates no 
variability; Pks < 0.005 indicates definite variability. Source XID 108 displays a multi-standard deviation flare. We extracted two spectra 
from the light curve of this source: one excluding events during the flare, and one including only flare events. 

















